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Abstract 

The weak phase 7 can be determined using untagged — > DKs or Bs — > D(f), Drj^'^ decays. In 
the past, the small lifetime difference y = Ar/(2r) has been neglected in B^, while the CP violating 
parameter e = 1 — has been neglected in both B^-B^ and Bs-Bg mixing. We estimate the 

effect of neglecting y and e. We find that in D decays to flavor states this introduces a systematic 
error, which is enhanced by a large ratio of Cabibbo-allowed to doubly Cabibbo-suppressed D 
decay amplitudes. 
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I. INTRODUCTION 



The interference between the decay chains B DK — > fj:,K and B DK foK is 
commonly utihzed in methods for extracting the phase 7 = arg (— V^dVjlVwjV^^) . The original 
idea, applied to these processes and to B^^-^ Dcf), was put forward in [l|, l2(]. Several variants 
of the idea can be found in , 0, 0, Slsl . The determination of 7 based on these methods 

is theoretically very clean. The uncertainty in the value of 7 is mainly due to low statistics 
of the relevant processes. Time-dependent decays B^ DKs studied in require flavor 
tagging of the initial -B", for which one pays a price in statistics. This is also the case for 



time-dependent Bg —>■ Dcj). Ref. discussed the extraction of 7 in untagged Bg decays, 
taking advantage of the finite nonzero width-difference between the two Bs mass-eigenstates. 
Recently it was shown that 7 can also be determined using untagged non-strange neutral 



B decays [12l], where the width-difference is too small to be measured. In this method, 
combining charged along with neutral B decays, and using several D decay modes, reduces 
the statistical error on the value of 7. 

The two parameters describing the width-difference and CP violation in mixing in the 

B^ systems. 



Ar ^ 



(1) 



are estimated to be below a level of 10~^ [l3|, [ij]. (Note that e is not a standard notation). 
Thus, usually these parameters can be safely neglected when dealing with untagged B^ — *• 
DKs- Indeed, both y and e were neglected in [1^]. 

In this work we consider the extraction of 7 from both untagged Bd DKs (here 
Bd = 5°) and untagged Bs Dcp, keeping explicitly y and e. We find that in D decays to 
flavor states, the sensitivity of 7 to y and e is enhanced by ratios of Cabibbo-allowed and 
doubly Cabibbo-suppressed D decay amplitudes. In such cases, neglecting y and e might 
result in a relatively large systematic error. This effect is expected to be more important in 
Bs decays, as the width-difference in the Bs system is expected to be much larger than in 



B. 
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15| . The effect in Bg decays may be taken into account through a measurement of the 
Bs width-difference. Eventually, when very large data sets will be available, a most precise 
determination of 7 will require using untagged Bd and Bs decays together with charged B 
decays. The purpose of this note is to point out the nontrivial enhancement of the effect of 
nonzero y and e, which may become relevant at that point. 

II. TWO-BODY D DECAYS 

Our analysis applies to both Bd — > DKs and Bs D(f), Drj^'^ . Working along the lines 



of 



121], we define for the first process, 

A{B' ^ D^Ks) = , A{B^ ^ D^Ks) = A„r„e*(^"+^) , (2) 
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where we take An and r„ to be real and positive, < (5„ < 27r and 7 are the strong and 
weak phases, respectively. Neglecting CP violation in K^-K^ mixing, amplitudes for CP- 
conjugate decays are given by inverting the sign of the weak phase 7. 

The parameter r„ is expected to be around 0.4 since it contains a Cabibbo-Kobayashi- 
Maskawa (CKM) factor |V^*bKs/KfeKis| ~ 0.4. It is not expected to be much smaller than 
this value unless color-suppression in 6 ^ ucs is more effective than in 6 — cus. An estimate 
of this effect can be obtained by comparing the ratios r* and r* measured in DK*^ 
and —>■ DK*^ decays. [The parameter r* is defined in analogy with ([2]), while r* is the 
corresponding parameter in charged B decays.] A recent BaBar analysis found an upper 
bound rl < 0.4 at 90% CL (fej, while a value r* = 0.11+^? was obtained by averaging a few 
measurements by BaBar and Belle 17|]. Using isospin and neglecting a pure annihilation 
amplitude describing B^ D^K*^, one has 121] 



, J r: = (4.0 ± 0.5)r: = 0.4+°-l . (3) 

Similar notations and considerations apply to Bs Dcj), Dri^'\ 

Considering D decays, we define amplitudes into a (possibly quasi) two-body final state 
fo and its CP-conjugate state /d by 

AiD^ ^ fo) = AiD' ^ fn) ^ Aj , A(D° ^ fo) = A{D^ ^ fo) = AfTje^'f . (4) 

By convention Af > 0, Vf > 0,0 < 6f < tt. We neglect CP violation related to D mesons, as 
well as D^-D^ mixing effects, which have been shown to be very small 18|] . The parameters 
rf,6f depend on the final state fo- While rf = 1 for CP-eigenstate (e.g. fo = K^K~) and 
Tf = 0{1) for non-CP fiavorless states (e.g. fn = K*^K~), Vf is very small for fiavor states 



(e.g. fo = K^TT^), Tf ~ 0.06 ~ tan^ [19ll- In this note we focus on a particular effect 
which applies to the latter class of final states. 

The amplitudes of the four decay chains are given by 

ABf = A(5° ^ foKs) = AnAf [rne'^'-+'''>rfe''f + l] , 
fnKs) = A^Aj[rfe''f +rne'^'--^^] , 
fnKs) = A^Aj[rfe''f +r^e'^'"+^^] , 
foKs) = A^Aj [r„e^('^"-^)r^e^^/ + l] , (5) 

with magnitudes 



M = 


A{B 


^Bf = 


A{B 


^Bf = 


A{B 




= A 


\M 


= A 


\^Bf 


= A 


\^Bf 


= A 



AnAf^ 1 + rlr) + 2r„rj cos {5f + 5„ + 7) 



n^fjrl + rf + 2r„r/cos(5/ - 5„ - 7) , 



AnAfJl + rlr) + 2r„rj cos (5/ + (5„ - 7) . (6) 



3 



We note that the phase of Asf/ABf depends also on 7. This dependence does not 
disappear in the hmit rj —>■ 0, in contrast to the dependence of the magnitudes of the 
amphtudes on 7. This will turn out to be a crucial point in our argument in the next 
section, when discussing the sensitivity of determining j to y and e for small values of rj. 

We now consider the time-dependent decay rate into fnKs for a state B^{t) which has 
evolved from an initially pure state. [A similar expression applies to T{Bs{t) — > 
Using our previously defined notations ([1]) and the standard notations 



T = Tt 



Am 



X 



r ' "'^pABf 

where a B flavor index [d or s) is implicit, one obtains [l3| . 



(7) 



r [B\t) ^ fnKs] 



\ ABf\^ [(coshyr + cosa;r) + |Af |^(coshyr — cosa;r) 
— 2(7?.e A/) sinhyr + 2(Xm A/) sinxr ] . 



(8) 



Similarly, for an initial B^ one has 



r [B\t) foKs] 



{Asff \p/q\'^ [(coshyr — cosxr) + |A/p(coshyr + cosxr) 



— 2{TZe A f) sinh yr — 2{Im A f) sin xr 
The untagged decay rate for a final foKs state is given by 
r [B\t) ^ foKs] + r [B\t) ^ foKs] = 



(9) 



— \ABf\^ { (1 + b/gp) [(1 + lA^n coshyr - 2(7^e A;) sinhyr] 



+ (1 — [(1 — |A/|^) cosxr + 2(Tm A/) sinxr] } 



(10) 



The total number of decays into foKs is obtained by integrating over time. To leading 
order in y and e, the time-integrated rate is 



— e 



Tf = \ABff + I^B/p - 2y \ABfABf \ cosipf 

X '^'^ 

\ABfABj\ sinv3/ 



X 



x^ 



2 (x2 + 1) 



A 



Bf\ 



A 



Bf\ 



(11) 



where we have defined = argA/. Similarly, for the CP-conjugate final state foKs one 
has 



/— \^Bf\ +1^-8/1 ~ 2y |^B/^_B/| cosy)/ 

2 



X 



X" 



\^Bf^Bf\ 



sm 



X 



2 (x2 + 1) 



A 



Bf\ 



A 



Bf\ 



(12) 
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Note that this result apphes regardless of whether the neutral B (Bg) mesons are produced 
hadronically or in e~^e~ collisions at T(45') [T(5S')]. In the latter case one neutral B [Bg) 
meson is observed decaying into foKs summing over all the decay modes of the 

second B [Bg). 

As expected, Eqs. (ITT!) and (fT2|) reproduce the result of [l2| in the limit y — 0, e — 0. 
We see that the entire x-dependence is suppressed by e. In the case of Bg decays, where 
a; ^ 1, we may expand in obtaining 

T/ = \ABf? + I^B/I^ - 2y \ABfABf\ COS^f - ^ (\ABff - I^B/T) 

+ O {y') + 0{e/x) + O (e^) . (13) 

Note that while the sum of the first two terms in Eqs. fill I) and (1121) depends on two 
combinations of B decay parameters, ^^(1 + r^) and cos (5„, the terms proportional to 
y and e involve An,rn and 5„ as three separate variables. This will affect the counting of 
parameters required for determining 7 which we discuss next. 

It was shown in |l3] that, when neglecting y and e, 7 can be determined from untagged 
B^ decays using several final states /|,. This holds true also when including nonzero values 
for y and e, where the number of unknowns is increased by one, as mentioned above. First, 
assume that y and e are known. Consider N different final states, {k = 1, . . . , N) 
which are not CP-eigenstates (e.g. = K~tt~^ , K*^K^), and their corresponding CP- 
conjugates, Assume that A^^^A^j and r'j-,r^ have been determined through branching 
ratio measurements in an independent sample of decays. For > 4, the A^ + 4 unknown 
parameters An, r„, 5„, 5^, and 7 can then be extracted by solving the 2A^ equations (JTTll and 
(fT2|) . (As a cross check for the solution one may use isospin in order to obtain estimates for 
An and r„ [JJ].) Including B° D*Ks followed by D* L)°7r°, D* L>°7, where the two 
strong phases 5* differ by vr 20], one has 6A^ measurements and A^ + 7 variables, a system 
which is solvable for N > 2. 

In principle, one may measure in this way not only 7 but also y and e. Including these two 
variables as unknowns permits determining 7, for A^ > 6 when using only the ground state 
D meson, and for N > 2 when using also decays involving D*. In practice, the values of y 
and e may be too small for a useful determination. The procedure above is expected to suffer 
from a large statistical error, and should therefore be used together with studying charged B 
decays. As we show in the next section, neglecting y and e simplifies the equations for neutral 
B decays, but introduces a systematic error which can be quite large for flavor-specific D 
decays. 



III. THE INDUCED THEORETICAL ERROR 

Let us now consider the error in 7 made by assuming y = e = 0, for the special case when 
D mesons decay to flavor states for which rj <^1. In order to estimate the sensitivity of 7 to 
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small nonzero values of y and e, we will take all other parameters as given (in the previous 
section we have outlined a procedure for obtaining a global solution for all parameters). 
Assuming that one measures the rate Tj for a certain decay mode, Eq. ffTTl) leads to a 
constraint of the form V^^'^ = Tf{'y; y, e). 

The dependence of F j on 7 in the limit y = e = is given by the first two terms on the 
right-hand-side of Eq. ( fTTjl . 

\ABff + \ABf\^ ocl + r2 +4r„rjcos5„cos(5/ + 7) + C(r5) , (14) 

where a proportionality factor A^^Aj was omitted. The effect of nonzero y and e is given by 
the remaining terms in Eq. (fTTj) . The terms in (fTT!) which are most sensitive to y and e are 
of zeroth order in rj, 

- 2y \ABfABf\ cosipf oc -2?/r„ cos((5„ - 2/5 - 7) + C(r/) , (15) 
ABfABfl siniff oc r„ sin((5„ - 2/5 - 7) + C(r/) . (16) 



Thus, while in f|T^ the term which depends on 7 is linear in r^, the 7-dependent terms in 
( ITSll and ( JT6II do not involve this small ratio. This implies that the sensitivity of determining 
7 to ?/ and e is enhanced by l/vf. The corrections to 7 from y ^ and e 7^ are given by 
terms of the form 

A7^ ~ ^ , AY ~ ■ (17) 



SO that the errors on 7 may be significant for fo = K^tt where r^+ji- — 0.06 19|. 



Using a recent evaluation of mixing parameters for the Bg system 15|, = 0.064 ±0.012, 



eg = (4.2 ± 1.2) ■ 10~^, and the measurement 2l| Xg = 26 ± 1, one notes that A7'^ can be 
safely neglected, while can be C(l) in Bg decays. Keeping the ?/-dependence when 
extracting 7 from time-integrated untagged Bg Dcj), Drj^'^ decays is thus mandatory. On 
the other hand, corrections in the B^ system, where one has measured [2l| x^ = 0.776±0.008, 



and one estimates [l5|] yd = (2.1 ±0.5) ■ 10~^, = (—9.6 ±2.2) ■ 10~^, are at a level of several 
percent. In this case one needs to consider corrections from y,e ^ when data reach high 
statistics. 

We stress that the above analysis provides only a rule of thumb for the size of the expected 
error in 7. The actual error in the analysis will depend crucially on the choice of the final 
states, and should be re-evaluated on a case by case basis. Large corrections from nonzero 
y and e apply only to flavor states in D decays, where r/ is small. Adding information from 
singly Cabibbo-suppressed D decays where rj ~ (9(1) is expected to dilute the error. 

IV. MULTI BODY D DECAYS 

Next we move to the case of multi-body D decays, such as D ^ Kst^^tt^. The error 
introduced by neglecting y and e may be enhanced by quasi two-body decays into flavor 
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states occurring in some parts of phase space, for instance in — > K*~^tt~. We first 
consider the model-independent approach presented in [sl, [l2| . We find that keeping y and e 
finite makes this approach non-practical, while the corrections can be implemented if the D 
decay amplitude is modeled by a sum of Breit-Wigner forms. We will focus on the three-body 
decay D Kstt'tt'^', however the following discussion is rather generic. 
Neglecting CP violation in D decays, one has 

A {D' Ks{pi),n-{p2),7r+{ps)) = A {D^ ^ Ks{pi),n-{p3),n+{p2)) 

= A(si3,si2)e*'(^i^'^-) , (18) 

where Sij = {pi +pj)'^. Denoting 

Aij = A{sii, sij), 6ij = 6{sii, sij) , (19) 
the amplitude of the cascade decay becomes 

Asf = A{B^ ^ foK) = A, [^326*'^^ + r„ e^*-^""^) A23 e^''^] , (20) 

and therefore 



ABf\= AnJAl^ + rlAl^ + 2v432A23r„ cos(7 + ^23 - + 5n) . (21) 



Similarly, 



AbA = AnJAl^ + rlAl^ + 2A23A32r„ cos(7 + ^23 - ^32 - 5„) . (22) 



These are the analogues of Eq. ([6]) . 

Using these magnitudes, the differential untagged decay rates, d?V j / dsi2 dsi^ and 
d^V j/dsi2dsi^, are given by Eq. (ITO!) . Dividing the Dalitz plot domain (si2,Si3) into bins, 
the time-integrated numbers of events over the ith bin are given by integrating Eq. (fTTi) over 
this bin. Thus, we have for the zth bin, 

N^ = Al {lit + I^Sr - 2yA^ + -^A|) + O [y') + O [e') , (23) 

where 

S+ = (1 + rfj7^+ + 4r„cos5n(ciCOS7 - Sisin7) , 
= (1 - r^)T^" - 4r„ sin 5„(ci sin 7 Si cos 7) , 

J i 

(24) 









/ dsi2 

J i 


dsi3 


AbjAb/I 


j dsi2 


dsis 


ABf^Bfl 
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and 



= [dsi2dsis {Al,±Al,) , 

J i 

Ci = jdSi2 dSi3 A23A32 COS (^23 " ^2) , 

Si = / dsi2 dsi3 A23A32 sin {S23 - S32) ■ (25) 

J i 

In order to see how many bins are needed if x, y and e were known, we consider m B- 
decay modes and divide the D-decay Dahtz plot into k bins. For each B decay mode one 
has three unknowns, y4„, r„, 5n- Assuming that the amphtudes are known from D decay 
data imphes four unknowns per bin, Cj, Sj, Af , vV?. If the bins are CP symmetric, the number 
of independent unknowns Cj, Sj reduces by half p], while no such symmetry reduction applies 
to A^, A|. Using both B^ DKs and B^ D*Ks, the number of unknowns including 7 is 
3k + 7, while the number of measurements is 3k. Such a system of equations is unsolvable. 
In principle, one could obtain sufficient information for the purpose of using known y and 
e to extract 7 from neutral B decays involving D Ks'n'^'n'^ alone. This would require 
including higher D resonances, or multi-body B decays [8|. However, this information and 
precise knowledge of y and e are not expected to be experimentally accessible very soon. 

We conclude that in the cases where one cannot neglect y and e (see the discussion 
in section [ITT]), the model-independent analysis of untagged multi-body -D-decays becomes 
impractical. A Breit-Wigner modeling of the density of events in the Dalitz plot j^, H, 23|. 
could be done if y and e were known. Without this a priori knowledge, the error in 7 could 
be large in regions of phase space where doubly Cabibbo-suppressed decays are important. 
A detailed study of the error would be needed. 



V. SUMMARY 

We have shown that the method presented in [l^ for extracting the CKM phase 7 
from untagged B^ decays leads to a systematic error caused by neglecting y and e. We 
found that the terms proportional to y and e are enhanced by a ratio of large to small D 
decay amplitudes. This enhancement is particularly notable when using doubly Cabibbo- 
suppressed D decays. It can lead to an error on 7 at the level of several percent in B^ —>■ DKs 
decays, and to an uncertainty of order one when extracting 7 in untagged B^ Dcp, Dr]'^'\ 

The method proposed in Ref. can thus be implemented without change in present 
B — > D^*^K^*^ data sets. However, future precision determinations of 7 will require taking 
into account corrections from y and e. While in Bs decays can be safely set to zero, the 
effect of a nonzero yg cannot be neglected, and experimental information of this parameter 
must be used already at the current level of precision. In the future, the discussed effects 
may have an impact on a global fit to 7 including both charged and neutral B decays. We 
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also note that at least in principle, with unlimited statistics, a measurement of y and e may 
be performed by using only time-independent experimental data. 
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